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ABSTRACT 

Context. The bright southern Cepheid RS Pup is surrounded by a circumstellar nebula reflecting the light from the central star. The 
propagation of the light variations from the Cepheid inside the dusty nebula creates spectacular light echoes that can be observed up 
to large distances from the star itself. This phenomenon is currently unique in this class of stars. 

Aims. For this relatively distant star, the trigonometric parallax available from Hipparcos has a low accuracy. A careful observation of 
the light echoes has the potential to provide a very accurate, geometric distance to RS Pup. 

Methods. We obtained a series of CCD images of RS Pup with the NTT/EMMI instrument covering the variation period of the star 
{P = 41.4 d), allowing us to track the progression of the light wavefronts over the nebular features surrounding the star. We measured 
precisely the phase lag of the photometric variation in several regions of the circumstellar nebula. 

Results. From our phase lag measurements, we derive a geometric distance of 1 992 ± 28 pc to RS Pup. This distance is affected by a 
total uncertainty of 1.4%, and corresponds to a parallax of ;r = 0.502 ± 0.007 mas and a distance modulus of /i = 1 1.50 ± 0.03. 
Conclusions. The geometric distance we derive is by far the most accurate to a Cepheid, and among the most accurate to any star. 
RS Pup appears both as somewhat neglected and particularly promising to investigate the mass-loss history of Cepheids. Thanks to its 
highly accurate distance, it is also bound to become an important luminosity fiducial for the long period part of the period-luminosity 
diagram. 

Key words. Stars: individual: RS Pup, Stars: circumstellar matter. Stars: distances. Stars: variables: Cepheids, ISM: reflection nebu- 
lae, ISM: dust, extinction 



^ 1 . Introduction - its pulsation period, which is one of the longest known values 

among Galactic Cepheids, is subjected to strong variations. 

5_j The variability of RS Pup was discovered by Miss Reitsma, 

C3 ■ Kapte yn's a ssistant in 1897, and independendy of her by I" spite of the wealth of information that can be retrieved 

■ " ' Innes (fT897] l. The period of about 41.3 days, determined at from thoroughly studying these peculiarities, RS Puppis has 

the turn of the 20* century, was the longest period among the been a neglected Cepheid for 2-3 decades. 

Cepheids known at that time. Long period Cepheids are im- Before our interferometric works (Kervella et al. l2006l 

portant because, owing to their larger mass, they evolve more Merand et al. 2006.. .2007J . the nebula around RS Pup was the 

rapidly than their shorter period counterparts. In addition, they only certain case of a direct relationship between a Cepheid 

are key objects in calibrating the period-luminosity relationship, and interstellar material. Another Cepheid, SUCas, also ap- 

and the exact knowledge of the behaviour of this relationship at pgars close t o a reflection nebula (Milone et al. I T9991 Turner & 

the long-period end is critical for establishing the cosmic dis- Evans lT984l l. but the physical association is more uncerta in. Th e 

tance scale because being more luminous, long-period Cepheids nebulosity around RS Pup was discovered by Westerlund ( IMTT i. 

can be discovered more easily in remote galaxies. No other Cepheid is known to be so intimately associated with 

RS Puppis (HD 68860, SAO 198944) is one of the most im- the interstellar matter. Westerlund already remarked that 'the 

portant Cepheids in our Galaxy due to its peculiarities: nebula may be large enough to permit the detection of possi- 

- it is embedded in a reflection nebula, ble variations in its intensity distribution due to the light varia- 

- it is situated in the region of a stellar association (Puppis III), tions of the star' . This offers a unique opportunity to determine 

the distance of the Cepheid by comparing the angular and lin- 

Send offprint requests to: P. Kervella ear radii of the denser rings in the reflection nebula. The po- 

* Based on observations collected at the European Southern sitional coincidence of RS Pup with the reflection nebula can 

Observatory, Chile under ESO Programme 078.D-0739(B). indeed be used for determining the distance and luminosity of 

Correspondence to: Pierre.Kervella@obspm.fr the Cepheid, based on the light echo phenomenon occuiTing 
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around such a large amplitude variable (Havlen [T972aJ see also 
Sugerman I2003I I. This distance is independent of the results 
obtained by other methods applicable for such variable stars. 
Another aspect of this positional coincidence is the study of the 
relationship between stellar pulsation and matter expelled from 
the star. 

Moreover, Westerlund (1963 ) pointed out that the reflection 
nebula surrounding RS Pup is a part of a stellar association con- 
sisting of early-type stars, which allowed him to derive the dis- 
tance of the newly detected Pup III association. The spatial co- 
incidence of the Pup III association and RS Pup allows another 
determination of the distance of the Cepheid independently of 
the methods based on the pulsation. 

Ring-like structures in the nebula surrounding the Cepheid 
imply that the circumstellar matter, at least partly, originates 
from the pulsating atmosphere of RS Pup. The mechanism for 
the formation of this circumstellar envelope (CSE) is currently 
unknown. Deasy (I1988J proposed several scenarios, ranging 
from evolution of the star before its Cepheid phase to mass 
loss during multiple crossings of the instability strip. As pointed 
out by Havlen (1972a) and Szabados ( |2003 |), it is unlikely that 
RS Pup is the only existing Cepheid-nebula association. Based 
on IRAS photometry and lUE ultraviolet spectra, Deasy ( 119881 ) 
identified mass loss in a number of Cepheids. The highest mass 
loss rate is attributed by this author to RSPup (10^^ Mg.yr"'). 
This level is very significant, comparable to that of the Mira stars 
(which host dense molecular and dusty envelopes), and explains 
the "bubble" structures that have been carved by this star in the 
interstellar medium. 

The available observational material on RS Pup mostly con- 
sists of photometric data covering more than a century. The 
analysis of these data is instrumental in determining the his- 
tory of variations in the pulsation period, which will be dis- 
cussed in a forthcoming paper. Changes in the pulsation period 
bear information on the evolutionary stage of the Cepheid (sec- 
ular period changes) and the physical processes occurring in the 
pulsating stellar atmosphere (period changes on time-scale of 
years). Recent high resolution spectroscopic observations with 
the HARPS instrument were reported by Nardetto et al. (120061 
120071 ). They show very broad spectral lines at certain phases, 
indicative of the presence of strong shock waves. 

In the present article, we describe our application of the light 
echo method to determine the distance to RS Pup. The observa- 
tions obtained with ESO's NTT are reported in Sect. |2] A de- 
tailed description of the processing steps that were applied to 
derive the photometry of the nebula is given in Sect. [3] Finally, 
we derive the distance to RS Pup in Sect. ID and compare our 
value to existing estimates in Sect. |5j 



2. Observations 

We observed RSPup using the 3.6m ESO New Technology 
Telescope (NTT) installed at La Silla Observatory (Chile), and 
equipped with the ESO Multi-Mode Instrument (EMMI). This 
instrument is a multipurpose imager and spectrograph (Dekker 
et al. II986I ). For the reported observations, we used the red 
arm of this instrument, based on two back-illuminated 2kx4k 
CCDs (15/im pixels) arranged in a mosaic, with a 7.82arcsec 
gap between the two detectors. The useful field of view is 
9.0x9.9 arcmin. We restricted our processing to 5.0x5.0 arcmin, 
still leaving a considerable margin around the nebula of RS Pup. 
In order to avoid a heavy saturation of the detector, we posi- 
tioned RS Pup in the gap between the detectors. The Nasmyth 



adapter of the telescope was rotated in order to align the direc- 
tion of the gap with a lower density part of the nebula. The pixel 
scale of EMMI is 0.1665 ± 0.0006".pixeri on the sk>Q, pro- 
viding a good sampling of the point spread function (PSF) even 
under good seeing conditions. 

Immediately before or after the RS Pup observations, we ob- 
served as a PSF reference the KlIII star HD 70555. We observed 
this star in addition to RS Pup in order to subtract the wings of 
the diffraction pattern of RS Pup from the nebula. We obtained 
a series of exposures through the B filter of EMMI at seven 
epochs, complemented by a single epoch of V, R, and jQ. The 
images in VRI colors will be discussed in a forthcoming arti- 
cle. The exposures were obtained randomly between 28 October 
2006 and 24 March 2007, but resulted in a reasonably good cov- 
erage of the pulsation cycle of the star The individual exposure 
times were relatively short, between 20 and 180 s, in order to 
avoid the saturation of the detector even very close to the star. 
The journal of the observations is presented in Table[T] The total 
shutter open time obtained on RS Pup with EMMI in all bands 
reaches more than 3.5 hours. The raw frames produced by the 
instrument were pre-processed in a standard way (bias removal, 
flat-fielding, cropping) using the Yoricli0 and IRAFQ software 
packages. The resulting frames were precisely co-aUgned using 
the IRAF script imcoadd (gemini package). 

In order to better vizualize the light echo phenomenon, we 
computed a continuous movie of the progression of the light 
echoes in the nebula, using a spline interpolation of the mea- 
sured epochs for each pixel (excepting the night of 20 December 
2006). Snapshots of this interpolated movie are presented in 
Fig. [1] for ten regularly spaced phases, and the movie itself 
(50 interpolated epochs) is available from the A&A web site in 
MPEG format. The variation of the luminosity of RS Pup is vis- 
ible on these images from the changing extension of the wings 
of its PSF. The morphological changes of the nebula with time 
are considerable, and they will be further studied in a forthcom- 
ing paper This movie was not used for the analysis developed 
below. 

To obtain absolute coordinates on a common world coordi- 
nate system (WCS), we retrieved from the IPAC archiv^ the 
FITS images covering the field of view of our images. The 
astrometric calibration was obtained by matching the coordi- 
nates of the brightest sources with their counterparts in the 
2MASS archive images. The absolute astrometric accuracy of 
the 2MASS images is about 0.1" (Ski'utskie et al. 12006] ). Due 
to the presence of stronger scattered light in the visible, our as- 
trometric uncertainty is estimated to be of the order of 0.2" (x; 
one pixel) over the EMMI field. As we do not aim at obtaining 
absolute astrometry, this calibration method is sufficient. 



' http://www.ls.eso.org/lasilla/sciops/ntt/emmi/newRed/ 
^ The transmission curves of the EMMI filters are available from 
!http://www.Is.eso.org/lasilla/sciops/ntt/emmi/emmLFiIters.html| 
"* http://yorick.sourceforge . net/ 

* IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

^ http://irsa.ipac.caltech.edu/ 
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Table 1. Log of the EMMI observations of RSPup and its PSF calibrator, HD 70555. MJD* is the median Modified Julian Date of 
the exposures. The phases (p were computed using the {Tq, P) parameters given in Sect. 13.21 



Date 


UTC 


MJD* 


Star 


Filter 


</> 


Exposures 


Airmass 


FWHM(") 


2006-10-28 


8:12:23 


54036.3419 


HD 70555 


B 




25s X 6 


1.17 


1.19 


2006-10-28 


8:37:14 


54036.3592 


RSPup 


B 


0.6854 


180 s X 10 


1.10 


1.14 


2006-10-29 


7:54:05 


54037.3292 


HD 70555 


B 




25s X 6 


1.21 


1.25 


2006-10-29 


8:23:17 


54037.3495 


RSPup 


B 


0.7093 


180 s X 10 


1.12 


1.13 


2006-12-20 


3:30:42 


54089.1463 


RSPup 


B 


0.9592 


180 s X 10 


1.40 


2.01 


2006-12-20 


3:57:31 


54089.1649 


HD 70555 


B 




25s X 6 


1.33 


1.97 


2007-01-13 


5:28:49 


54113.2283 


RS Pup 


B 


0.5404 


180 s X 10 


1.01 


0.91 


2007-01-13 


5:54:39 


54113.2463 


HD 70555 


B 




25s X 6 


1.01 


0.70 


2007-01-22 


2:39:17 


54122.1106 


HD 70555 


B 




25s X 7 


1.16 


1.16 


2007-01-22 


3:09:38 


54122.1317 


RS Pup 


B 


0.7552 


180 s X 12 


1.08 


0.96 


2007-01-27 


3:49:50 


54127.1596 


HD 70555 


B 




25s X 6 


1.02 


0.88 


2007-01-27 


3:55:38 


54127.1636 


HD 70555 


V 




6s X 6 


1.02 


1.13 


2007-01-27 


4:00:11 


54127.1668 


HD 70555 


R 




2s X 6 


1.01 


0.98 


2007-01-27 


4:04:58 


54127.1701 


HD 70555 


I 




3s X 6 


1.01 


0.98 


2007-01-27 


4:18:53 


54127.1798 


RSPup 


B 


0.8770 


179s X 5 


1.01 


0.88 


2007-01-27 


4:35:02 


54127.1910 


RS Pup 


V 


0.8773 


45s X 8 


1.00 


0.89 


2007-01-27 


4:47:06 


54127.1994 


RSPup 


R 


0.8775 


20 s X 12 


1.01 


0.80 


2007-01-27 


4:59:11 


54127.2078 


RSPup 


1 


0.8777 


30 s X 10 


1.01 


0.81 


2007-03-24 


0:16:25 


54183.0114 


RSPup 


B 


0.2244 


180 s X 10 


1.01 


0.88 


2007-03-24 


0:44:32 


54183.0309 


HD 70555 


B 




25s X 6 


1.00 


0.74 



3. Data analysis 

3.1. PSF subtraction and photometry 

3.1.1. Preparation of the reference PSF 

On the images of RS Pup, large spikes are created by the sec- 
ondary mirror spider support and the inhomogeneities of the op- 
tical surfaces in the telescope. Fortunately, from one star to the 
other and over long periods of time, the distribution of light in 
this diffraction pattern is very much reproducible, and it can be 
efficiently measured and subtracted. We therefore observed im- 
mediately before or after RS Pup the bright KlIII star HD 70555. 
As shown in Fig. |2] (top), the field around this star is relatively 
rich in background stars, due to its low Galactic latitude (+2°). 
Its position near RS Pup is an advantage to match the airmasses, 
but the presence of background stars, which must be removed 
before subtraction from the RS Pup images, creates some diffi- 
culties. 

Due to the alt-az configuration of the NTT, the pupil of the 
telescope (and all the attached PSF defects) rotate at an angular 
rate that is a function of the altitude and azimuth of the object. 
The position angle of the diffraction spikes is given by the paral- 
lactic angle at the time of the exposure. These optical defects are 
caused by the secondary mirror support and the optical elements, 
and therefore they rotate with the telescope pupil. During nor- 
mal operation of the telescope, the field de-rotator (a mechanical 
device on which the instrument is mounted) takes care of rotat- 
ing continuously the instrument in order to keep the position of 
the objects constant on the detector. Of course, any defect intro- 
duced by the instrument itself does not rotate with the telescope 
pupil (such as ghosts created by internal reflections, e.g. inside 
filters). 

To take advantage of this effect, we "re-rotated" the PSF 
exposures obtained on HD 70555 for the nights of 2006-10-28, 
2007-01-13, 2007-01-27 and 2007-03-24. In other words, using 
the parallactic angle information, we aligned precisely the PSF 
spikes with each other, letting the background star field rotate. 
We chose these four nights as their seeing FWHM is compa- 
rable (from 0.70" to 1.14"), and they are all different in terms 
of parallactic angle. The effective FWHMs of these exposures 



were then all equalized to 1.19", using the psfmatch task of 
IRAF. This command convolves the images by a kernel that is 
computed to match precisely the images to a reference frame 
(in this case, the first observation of HD 70555 obtained under 
slightly degraded seeing conditions on the night of 2006-10-28). 
Computing the median of these re-rotated images resulted in a 
very efficient removal of the background objects, leaving only 
the PSF of the central star, as shown in Fig. |2] (bottom). The 
latter image is rotated to a null parallactic angle, and the im- 
age scale is logarithmic to show the faint extensions of the PSF 
wings. The background stars disappeared almost completely. As 
in all the images presented in this article, the position of the ob- 
served star is marked with a spot. Apart from the small artefacts 
close to the center, this image shows the actual PSF produced by 
the NTT with a seeing of 1.19", for a null parallactic angle, i.e., 
with the altitude direction of the telescope aligned along the ver- 
tical axis of the image. We finally prepared an adapted reference 
PSF frame for each individual image of RS Pup by re-rotating 
the image presented in Fig. |2] (bottom) to the same parallactic 
angle. 

3.1 .2. Preparation of tine images of RS Pup 

In order to properly subtract the image obtained at Sect. 13.1.11 
we matched the PSFs of each RS Pup image to the same FWHM 
of 1.19", which is also the FWHM of the PSF reference im- 
age. We normalized photometrically all images to the first one 
using a sample of twenty field stars and aperture photometry 
(with a large 20 pixels aperture) with the phot task of IRAF. 
The measured standard deviation of the photometry obtained on 
these twenty stars using a 5 pixels aperture is ±0.03 mag over 
our complete image sample. 

The night of 20 December 2006 is a special case as the very 
large FWHM is not caused by a poor seeing, but by a defocus 
of the telescope (the true seeing was around 0.7"). In order to 
recover useful information from this night, we used a special 
processing method. For this particular night, we degraded the 
FWHM of the reference PSF image from Sect. |3.1.T] to match the 
2.01" FWHM of the RS Pup image. We thus obtained a special 
PSF image set matching this series of images. 
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1 arcmin 



Fig. 2. Top: Image of the field around HD 70555 showing the 
wings of the PSF of the NTT and the background stars. Bottom: 
Reference NTT/EMMI PSF extracted from the observations of 
HD 70555 obtained over four epochs. 

3.2. PSF subtraction and photometry of RS Pup 

The goal of the whole procedure described in this section is to 
subtract the wings of the RS Pup PSF without affecting the neb- 
ular features over which the light echoes propagate. For this pur- 
pose, it is necessary to compute the flux ratio between RS Pup 
and HD 70555. We initially thought about computing this ratio 
from predicted magnitudes of RS Pup, but in practice this is not 
possible due to its irregularly evolving period. This makes it dif- 
ficult to phase precisely our measurements over archive photom- 
etry (even if only a few years old). So we resorted to a direct 




* 
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Fig. 3. Top: Rectangular zones (shown as light grey patches) se- 
lected to compute the photometric ratio between RS Pup and the 
reference PSF frame, superimposed on a single frame obtained 
on 28 October 2006. Bottom: Result of the PSF subtraction of 
the top frame. 



measurement of the flux in the wings of the PSF of RS Pup. 
However, the measurement itself is made difficult by the pres- 
ence of the relatively bright nebular features that tend to bias the 
photometry towards higher flux. 

We thus selected two regions close to RS Pup that show a low 
nebular contribution on all images (Fig. [3] top). Over these two 
regions, we smoothed the result of the division of each image 
of RS Pup by the corresponding PSF image, using a median box 
filter We then took the minimum value as the ratio between the 
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Table 2. Measured photometry of RS Pup. 



Date 


MJD 


lJIld.sC 


fllg 




2006-10-28 


54036.3592 


0.6854 


9.329 


0.033 


2006-10-29 


54037.3495 


0.7093 


9.343 


0.038 


2006-12-20 


54089.1463 


0.9592 


1.11% 


0.059 


2007-01-13 


54113.2283 


0.5404 


8.970 


0.040 


2007-01-22 


54122.1317 


0.7552 


9.315 


0.044 


2007-01-27 


54127.1798 


0.8770 


8.564 


0.037 


2007-03-24 


54183.0114 


0.2244 


8.328 


0.067 



two frames. This method ensures that we do not overestimate the 
flux from RS Pup due to the nebula. As an example, the result of 
the subtraction is shown in Fig.[3](bottom) for a single frame. For 
each night, we computed the median of all subtracted images, an 
operation that further reduced the residuals of the PSF subtrac- 
tion. As an example of the result of this processing, the image 
of RS Pup obtained on 27 January 2007 is shown in Fig. |4] As 
visible in this image, the field de-rotator of the NTT was off'set 
by 60° for the observations to avoid masking dense parts of the 
nebula inside the CCD gap. North is up and east is to the left on 
this image, but the other figures of this article are not corrected 
for the rotational ofi'set. 




Fig. 4. PSF-subtracted image of the RS Pup nebula obtained on 
flie night of 27 January 2007 at a phase </> = 0.8770. 



The derived flux ratio is the same between the frames ob- 
tained over one night to a few percent. We used the standard de- 
viation of this dispersion as the statistical uncertainty of the ratio. 
To estimate the systematic error, we obtained aperture photom- 
etry of a sample of faint stars close to RS Pup (but unaffected 
by the nebula), in the subtracted images. The dispersion of the 
measured magnitudes at all epochs (+0.03 mag) was taken as the 
measurement uncertainty. 

The PSF calibrator (HD 70555) was chosen in the spec- 
trophotometric catalogue assembled by Cohen et al. ( 1999). This 
gives us the possibility to use it as a convenient photometric ref- 



erence to convert the measured ratio with RS Pup into magni- 
tude. The Hipparcos catalogue (ESA'1997) gives niB - 6.249. 
Although neither star was actually present on the detector, the 
measured photometry of RS Pup is relatively accurate (Table |2] 
and Fig.|5]). From these magnitudes, we derive the following ref- 
erence epoch of the maximum light: To = MJD 54090.836+0.08 
(JD 2454091.336). For this fit, we adjusted a phase shift A0 to 
match the measured photometric data points, using an interpo- 
lated light curve constructed from archival photometry (using 
Berdnikov's dat^. We find no systematic magnitude shift down 
to a ±0.02 mag level. 

While the shape of the light curve has remained stable 
(which justifies the fit of the recent photometric data to the nor- 
mal light curve), the period of RS Pup has varied dramatically 
over the last century. There are short term changes superimposed 
on a uniform secular period increase (to be discussed in a forth- 
coming paper). The instantaneous period of f = 41.4389 days 
was derived from the most recent part of the "Observed minus 
Calculated" {O - C) diagram, i.e., from the O - C residuals cal- 
culated for the epochs between 2002 and 2006. We used this 
(To, P) combination to phase our measurements. 



7,5 




9,5 \ . . . . ^ 

0,0 0,2 0,4 0,6 0,8 1,0 

Phase 



Fig. 5. Photometry of RS Pup from our observations, superim- 
posed on its B band hght curve for - MJD 54090.836 and 
P = 41.4389 days. 



After the PSF subtraction, the images of the night of the 20 
December 2006 were deconvolved using the Lucy-Richardson 
algorithm in order to match the PSF FWHM with the other im- 
ages in the series (normalized to a FWHM of 1.19"). After this 
procedure, we checked the consistency of the photometry of a 
sample of stars in the field and found a good agreement within 
less than +0.1 mag. Although the quality of the data of this par- 
ticular night is low, this step is important in the context of neb- 
ular photometry to limit the contamination among neighbouring 
parts of the nebula. 

The variation of the aspect of the nebula over the six phases 
for which we have high quality images (excepting the night of 
20 December 2006) is presented in Fig.|6] 



available from |http://www.sai.msu.su/groups/cluster/CEP/PHE/| 
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4. Distance determination 

4.1. Selection of the nebular features 

The determination of the distance to RS Pup following the 
method of Havlen ( |1972a| ) is based on the measurement of the 
phase difference between the variation of the star and the varia- 
tion of isolated nebular features. When a particularly dense and 
localized feature can be isolated, we can make the reasonable 
hypothesis that the light diffused by this feature dominates the 
light diffused by the other superimposed layers. In other words, 
although several layers contribute to the local photometry (at the 
same projected position on the plane of the sky), the phase offset 
of the photometric variation can be considered identical to that 
of the localized "blob". 
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Fig. 7. Positions of the isolated nebular features chosen in the cir- 
cumstellar nebula of RS Pup. Due to their phase-shifted variation 
cycles, not all features are simultaneously visible on this image, 
obtained on 13 January 2007 at a phase of = 0.5404. The size 
of the circles marking each region corresponds approximately to 
the 5x5 pixel boxes used to average the derived parameters. 

The selection of the regions of interest was done in two steps: 

1. identify visually all the well-defined compact nebular 
features that show a clear variation over several images, 

2. select within this sample the features showing the largest 
photometric variation amplitude. 

Step 2 is necessary to select the features that are most hkely to 
correspond to a single position in the nebula, thus avoiding the 
confusion of several superimposed layers located at different dis- 
tances. This also selects the most accurate phase measurements. 
This procedure resulted in a selection of ten prominent nebular 
features. Figure|7]shows the position of these regions. We chose 
features that are located as close as possible to the star, although 
avoiding the closest regions where residuals from the PSF sub- 
traction are present. They are the most discriminating for the dis- 
tance determination, as the integer number of cycles A^, between 



their variation curve and the curve of the Cepheid is small. This 
translates into a large shift in distance for different values of A^, . 

4.2. Photometry of the nebular features 

The measurement of the surface magnitude of an inhomoge- 
neous nebula depends on the choice of aperture size. Choosing 
a too large aperture can lead to a degraded estimate of the vari- 
ation of each individual component, while a too small aperture 
will produce inaccurate measurements due to the limited photo- 
metric signal-to-noise ratio (SNR). One advantage of the careful 
normalization of the PSF of our images before their photometric 
co-referencing is that our nebular photometry is mostly insensi- 
tive to the seeing fluctuations between each night. We can there- 
fore use each pixel of the nebula image individually to measure 
photometry, when the SNR is sufficient. We thus did not adjust 
centroids on the nebular features, as they are highly variable, 
but instead we chose one fixed position for each of them over 
which we obtained the photometric measurements. The effective 
photometric aperture on each pixel therefore corresponds to the 
seeing PSF, with a constant FWHM of 1.19" for all images (see 
Sect. [Trill- 
in all the subsequent analysis of the nebula variation, the 
night of 20 December 2006 will be used only for verification 
purposes, and not for the fit itself. The deconvolution procedure 
that we used specifically for this night shows no indication of 
introducing a photometric bias on the photometry of the stellar 
sources, but could possibly create artefacts in the surface pho- 
tometry of the nebula. 

4.3. Measured phase offsets 

The measured properties of the selected features and the am- 
plitude and phase of their respective photometric variations are 
listed in Table[3] 

The phases and amplitudes were subsequently averaged over 
a 5 X 5 pixels box (0.83" x 0.83") centered on each position. 
The fit was obtained by adjusting the phase, mean value and 
amplitude of the photometric curve of RS Pup to the photom- 
etry of each pixel using a classical minimization through a 
Levenberg-Marquardt algorithm. The fit was obtained separately 
on each pixel of the 5x5 pixel boxes. The standard deviation of 
the measurements over the 5x5 pixels box, added quadratically 
to the fitting error bars was taken as its associated uncertainty. 
While the estimation of each phase offset is more precise, this 
conservative approach is intended to account for possible photo- 
metric contaminations by other nearby features. The final listed 
phase and amplitude values correspond to the average of these 
parameters over these boxes. Figure [8] shows the best fit curves 
for the selected nebular features detailed in Table [3] In this fig- 
ure, the plotted error bars correspond to the standard deviation of 
the 5x5 pixel window used to average the adjusted parameters. 
The data point marked with a cross was excluded from the fit- 
ting process as it corresponds to the night of 20 December 2006, 
when the telescope was accidentally out of focus. 

It is interesting to remark that the feature #1 selected by 
Havlen ( |1972a| l corresponds approximately to our feature #2, al- 
though an exact identification is made difficult by the absence 
of relative coordinates in his article (we used his Fig. 1). Havlen 
obtained a minimum light phase of (p - 0.52 - 0.57, depend- 
ing on the level of dilution of the photometric curve, while we 
obtain ^ 0.55 (Fig. |8]i. Considering the unknown uncertainty 
on Havlen's feature position, the agreement is satisfactory. The 
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Table 3. Measured photometry of the selected nebular features. The x and y columns list the coordinates of the selected pixels in 
the data cube. The angular radius 6j of these pixels is counted in arcseconds from RS Pup, and the azimuth is measured from north 
(A^ = 0°, £■ = 90°). The coordinates of RS Pup in the image, in pixels from the bottom left, are x = 619.3 + 0.2, y = 628.2 + 0.2. 
"Mean," is the average value, and A, the amplitude of the adjusted photometric curve, expressed in ADU and A0, is the phase 
offset with respect to RS Pup, with its associated total error (see text). A^, is the computed integer number of periods between the 
variation cycle of RS Pup and that of the feature, c/, is the derived distance to RS Pup. The six selected features for the final distance 
computation are marked with * in the first column (see Sect. |4.4| for details). 



# X (pix) y (pix) 9i {") Azimuth (°) Mean, (ADU) A, (ADU) A^, Nj dj ± a (pc) 



1 


578 


748 


21.10 


139 


185 ±3 


313 ± 11 


0.983 


± 0.020 


5 


2034 ± 


7 


2* 


522 


710 


21.16 


170 


185 + 15 


280 ± 17 


0.809 


±0.012 


5 


1969 ± 


5 


3* 


493 


659 


21.65 


196 


175 + 19 


274 ± 43 


0.989 


±0.013 


5 


1985 ± 


5 


4* 


528 


668 


16.58 


186 


171 ± 11 


241 ± 19 


0.576 


± 0.009 


4 


1980 ± 


4 


5* 


549 


694 


16.03 


167 


232 ± 14 


275 ± 28 


0.504 


± 0.083 


4 


2016 ± 


38 


6 


542 


624 


12.89 


213 


181 ±3 


291 ± 20 


0.102 


± 0.023 


3 


1726 ± 


14 


7* 


568 


587 


10.96 


249 


155 ±9 


274 ± 15 


0.098 


± 0.023 


3 


2028 ± 


16 


8* 


655 


456 


29.28 


312 


99 ±9 


151 ± 20 


0.207 


± 0.036 


8 


2010 ± 


9 


9 


695 


717 


19.42 


80 


87 ±6 


83 ±4 


0.697 


±0.016 


5 


2104 ± 


6 


10 


569 


719 


17.28 


149 


197 ± 6 


220 ±9 


0.602 


±0.014 


4 


1910 ± 


6 



three other features he selected are located further away from 
RS Pup and do not correspond to ours. 

4.4. Distance estimation and statistical uncertainty 

We assume in the following that the integer number of periods 
A^, separating the photometric variation of RS Pup from the vari- 
ation of each nebular feature is increasing with the angular dis- 
tance from the star. This is a natural consequence of our assump- 
tion that the selected features are located close to the plane of 
the sky (see discussion above). This means that a nebular feature 
with index / located closer to RS Pup than another region j will 
have A^, < Nj. From the phase offset A0, of each feature (with 
1 < ! < 8), we can derive a distance J, to RS Pup assuming an 
arbitrary value for A^,, from the simple relation: 



(Ni + Acf>i)P 
5.7755 10-3 0: 



where d, is in parsecs, P in days, and 0, in arcseconds. As rightly 
pointed out by Havlen ( |1972al l, it is not possible to obtain a 
unique distance solution, independently of the number of mea- 
sured features. However, for the right values of the A^, numbers, 
we expect a marked convergence of the estimated c/, towards a 
single value. As we assumed that A^,- is increasing with the sepa- 
ration from RS Pup, we can summarize the determination of the 
distance to a study of the relative standard deviation of the d/ dis- 
tances as a function of the a priori assumed d distance to RS Pup. 
Figure |9] clearly shows the best convergence for a starting value 
around 1 990 pc. The best-fit unweighted average distance d of 
all di in this case is d - I 992 + 20 pc. If we apply a weight- 
ing of the dj measurements by the inverse of the phase variance 
of the corresponding nebular feature, we obtain 1 983 pc, within 
half a cr of the previous value. The unweighted median distance 
is 1997pc. 

The individual formal error bar on each di value is usually 
below 1%, and in some cases even better than 0.3%, but their 
dispersion is larger due to their non-planar spatial distribution 
around RS Pup. To take this into account, the uncertainty on 
the average distance was computed using the unweighted boot- 
strapped standard deviation of all di, minus the two highest and 
two lowest values of our set. A description of the bootstrapping 
technique can be found in Appendix B of Kervella et al. (2004). 
This clipping of four di distance values (out of ten) is useful to 



reduce our sensitivity to the dispersion of the measured blobs rel- 
atively to the plane of the sky. Features located closer or further 
from us than RS Pup will have a smaller apparent angular sepa- 
ration from RS Pup, and the distance estimate will be uncertain. 
We suspect this is the case for the nebular feature #6 in Table [3] 
However, fits with all ten blobs included, or inversely with six 
features clipped out, leads to exactly the same unweighted aver- 
age value of 1 992 + 19 pc, so this clipping does not introduce 
any significant bias. 




500 1000 1500 2000 2500 3000 

A priori starting distance (pc) 

Fig. 9. Relative dispersion of the di distance estimates as a func- 
tion of the initial assumed distance to RS Pup for 10 (thin black 
curve) or 21 (thick grey curve) nebular features. The arrow 
points at the minimum dispersion of 1 . 1 %. 



As shown in Fig. |9] other starting distances than 1 990 pc 
provide some convergence, in particular around 650 pc and 
1 340 pc, but the dispersion of the di is much larger than for 
1 992 pc (2.4% and 2.8%, compared to 1.1%). They correspond 
to different sets of increasing A^,- values. The corresponding re- 
duced are respectively - 73, 775 and 151, for the three 
increasing distances. The smaller for 650 pc is due to two 
effects: 
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10 15 20 25 

Radial distance from RS Pup (arcsec) 

Fig. 10. Total phase A^, + 0, of the observed variations of the neb- 
ular features around RS Pup as a function of their radial angular 
distance from the star for the A^,- values given in Table [3] The 
line is a two-parameter linear fit through the points (see text). 
The data points shown as open circles were excluded from the 
fit (two largest and two smallest distance estimates from the set). 



4.5. Systematic error sources and final distance 
4.5.1 . Astrometry and photometry 

The pixel scale of the EMMI detector is 0.1665 + 0.0006".pix ^ 
resulting in a systematic uncertainty of 0.4% on all 0i values, 
and consequently on the distance. The position of RS Pup was 
determined with an accuracy of +0.2 pixel on the reference im- 
age, from the intersection of the large diffraction spikes created 
by the secondary mirror support (visible in Fig.|2] bottom). This 
translates into a relative uncertainty of 0,033" on each measured 
Oi. For the smallest angular distance of feature #7 (0, = 10.96"), 
this corresponds to a relative uncertainty of 0.3%. 

The angular separations 0, themselves do not contribute to a 
systematic uncertainty as we considered pixel to pixel photome- 
try, without centroiding the features. The accuracy of the regis- 
tration of the images with respect to each other is very high as 
we used a cross-correlation over hundreds of background stars 
(star positions are constant on the detector). The photometric 
normalization of each image with respect to the first in the se- 
ries was achieved using a sample of field stars with an accuracy 
of a few 0.01 mag (Sect. [3TL2l i. The subtraction of the telescope 
PSF introduces additional photometric errors, but they are taken 
into account in this statistical uncertainty. More generally, the 
photometry is not expected to be a significant contributor to the 
overall systematic uncertainty, as the residual fluctuations are re- 
flected into the statistical phase uncertainties listed in Table [3j 



- The error bars on each phase offset translate into constant 
uncertainties on the di values, independently of the chosen 
A^, . If the average distance becomes smaller (by reducing all 
the A^, simultaneously for instance), the relative error of each 
di increases, while the errors themselves remain constant. As 
a consequence, for the smaller average distance, the be- 
comes smaller, while the standard deviation of the c/, remains 
the same. The is therefore a poor indicator of the con- 
vergence of the method, compared to the straight standard 
deviation. 

- The 650 pc distance is nearly a third of 1 992 pc, and it be- 
trays a geometrical aliasing effect of the distribution of the 
selected blobs on the nebula. We can exclude this small dis- 
tance (as well as 1340pc) using our verification fit with 
a distinct set of 21 nebular features (see also Sect. I4.5.4l i. 
which gives x^^^ - 924 around 650-700 pc compared to 
^id = 674 around 2 000 pc. 

One should also note that a distance of exactly twice 1 992 pc 
is also a possibility. This is a natural consequence of the method, 
as the number A^, of cycles is exactly doubled for each blob, but 
such a value of 4000pc is excluded by previous determina- 
tions from other methods (see Sect.jSj. This also applies to larger 
integer multiples. 

We can check a posteriori that the derived distance is self- 
consistent by fitting the derived total phases A^, + 0, as a function 
of the angular radius from RS Pup for both the slope and the zero 
point. The slope is proportional to the inverse of the distance, and 
the zero point should give a null phase. The derived distance is 
1 993 pc, and the residual phase is very small {A<po - 0.007) for 
the selected features (Fig.fTOli. The agreement is not so good for 
the other distances where a local convergence is observed, giving 
further confidence that they are not correct. 



4.5.2. Period and reference epoch 

The computation of the phase offset relies on the assumed (P, Tq) 
pair An error on these two values will have an impact on the 
determined phase offset of a given feature, and thus on its esti- 
mated linear distance from RS Pup. We derived the time of the 
actual maximum light epoch Tq from our own data, and thus 
obtained a value suitable for the epoch of our observations. An 
uncertainty of 0.08 day is attached to this estimate. It translates 
into a maximum systematic uncertainty of 0.08/41.4 = 0.2% on 
the distance. As they are spread over 147 days, our observations 
cover about 3.5 periods of RS Pup. Assuming an uncertainty of 
+0.02 day on the period P, this will create a maximum ±0.07 day 
uncertainty on the last maximum dates, equivalent to +0.2% on 
the distance. 

4.5.3. Integer number of periods A^, and scattering 

The computation of the distance for each nebular region is based 
on the hypothesis that they are located in the plane of the sky, 
at the same distance as RS Pup itself. This is underpinned by 
the reasonable assumption that the distribution of the features 
around the star is, on an average, isotropic. Several features 
could be located closer or farther from us than RS Pup, result- 
ing in a larger or smaller derived distance from these particular 
regions, depending on the integer rounding error made on the A^, 
values. However, it is unlikely that it is simultaneously the case 
for all features, especially as we consider features located close 
to the star, hence with low A^, values (Table O. Considering that 
this uncertainty due to A^, is of statistical nature (averaging out 
with several blobs), we do not introduce here a specific system- 
atic error contribution. The statistical dispersion of the measure- 
ments is included in the bootstrapped error bar given in Sect. 14.41 
The scattering angle of the light from the Cepheid by the dust 
can in principle create an uncertainty on light echo distances. 
This question is discussed extensively by Havlen ( |1972bl l. 
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Although this is true in the case of scattering by a continuous 
dusty medium, we rely in the present work on spatially isolated, 
well defined nebular features. This makes our distance estimate 
mostly insensitive to the scattering law. At the limit, the scat- 
tering by an unresolved, isolated compact dust region would be 
fully independent of the scattering function. For this reason, we 
do not introduce a specific systematic uncertainty due to the un- 
known scattering function. 

4.5.4. Choice of the nebular features 

The choice of the nebular features and their number could pos- 
sibly create a bias. As a verification, we considered a distinct 
sample of 21 positions in the nebula, extending up to 0, - 40" 
from RS Pup. We derived from this sample an unweighted aver- 
age distance of2025+ 18pc (Fig. |9] thick grey curve), only Icr 
away from the distance obtained in Sect. |4.4| (l 992 + 20pc). 

The slightly different distance value is due to the inclu- 
sion of poorly defined nebular features that show a "smoothed" 
photometric variation curve due to the superimposition of sev- 
eral layers located on the line of sight. The estimation of the 
phase for these features is thus slightly biased, as shown by 
Havlen (T972ai An indication that this bias is present is that the 
unweighted median distance from this large sample is 2 012 + 
22 pc, closer to the value obtained with the selected features. 
Moreover, the zero point of the retrospective two-parameter fit 
of the blob distances gives - -0.10, significantly worse than 
for our carefully selected reference sample. 

From this result, we affect a systematic uncertainty of +0.8% 
(±16pc) to the choice of the nebular features. 

4.5.5. Final distance and uncertainty 

As discussed in the previous section, several systematic uncer- 
tainty sources have an impact on the distance of RS Pup deter- 
mined from its light echoes. Table |4] lists the most important 
ones, with their estimated relative impact on the determined dis- 
tance. All these terms have a linear effect on the measured dis- 
tance. Adding quadratically all these contributions, we obtain a 
total systematic uncertainty of 1.0% (Table |4|. This small num- 
ber shows that the light echo method is very robust in the case 
of a periodic star with a well-defined variation curve. In partic- 
ular, the selection of the nebular features, although an important 
point, is not a strong limitation to the final accuracy. Writing 
separately the statistical and systematic uncertainties, the final 
distance to RS Pup is therefore: 

d = 1992 + 20 + 20pc. 

corresponding to a parallax of tt = 0.502 + 0.007 mas (total 
relative uncertainty of 1.4%) and a distance modulus of // = 
11.50 ±0.03 mag. 

5. Comparison with existing distance estimates 

5.1. Trigonometric parallax 

Due to its large distance, the trigonometric parallax of RS Pup 
measured by Hipparcos cannot be accurate. The original 
Hipparcos catalogue published a value of tt = 0.49 + 0.68 mas 
(ESA 1 19971 1, which is not significant. However, a revised par- 
allax was published by van Leeuwen et al. (2007) as tt = 
1.44 + 0.51 mas, corresponding to a 3cr detection, and the fi- 
nal value that appears in the revised Hipparcos parallax cata- 
logue is TT = 1.92 + 0.65 mas (van Leeuwen i2007ai see also 



Table 4. Systematic uncertainty contributors on the distance of 
RS Pup. 









Pixel scale 


0.1665 ±0.0006" 


0.4% 


Position of RS Pup 


±0.2 pixel 


0.3% 


Reference epoch 


MJD 54090.836 ±0.08 


0.2% 


Photometric period 


41.4389 ±0.02 days 


0.2% 


Choice of features 


±16pc (see text) 


0.8% 


Total 


±20 pc 


1.0% 



van Leeuwen I2007bl l . Given the limited accuracy of Hipparcos 
(~ 1 mas), we do not have an explanation for this apparent de- 
tection of parallax, in disagreement with our determination. 

5.2. Infrared surface brightness 

A more promising distance determination comes from the 
infrared surface brightness method, a variant of the Baade- 
Wesselink technique. Barnes et al. (2005) included RSPup in 
their catalogue, with a Bayesian distance of 1 984 + 122 pc 
and a least squares distance of 2 004 ± 59 pc, both in excellent 
agreement with the present determination, although less accu- 
rate. However, a systematic uncertainty arises with this tech- 
nique (not accounted for in the above error bars), due to the 
projection factor used to convert radial velocities to pulsation 
velocities. Barnes et al. used a classical p-factor first derived 
by Gieren et al. (1993), namely p = 1.39 - 0.03 logP, which 
amounts to 1.341 in the case ofRS Pup (log P = 1.6174). Fouque 
et al. (I2007I I derived a distance of 1 830 pc based on the same 
data, but using a lower projection factor based on a mean rela- 
tion derived by Nardetto et al. ( |2007l ), p = 1.366 - 0.075 logP 
(cross-correlation method), which leads to 1 .245 in the case of 
RS Pup. This revision is based on recent hydrodynamical models 
of RS Pup itself, combined with high-resolution spectroscopic 
measurements. The question of the /9-factor will be further dis- 
cussed in a forthcoming paper. The interested reader can refer to 
Groenewegen's (120071 ) recent work. 

In conclusion, we can regard the infrared surface brightness 
distance as supporting our result, although both its random and 
systematic uncertainties are larger than those of the present ge- 
ometric distance measurement. For the sake of better visualiza- 
tion. Table |5] gives a summary of the distance values of RS Pup 
obtained from various observational methods. 

5.3. Stellar association 

RS Puppis was found to be a likely member of the Pup III OB as- 
sociation discovered by Westerlund (1963). Based on the color- 
magnitude diagram of about two dozen early-type stars in the 
region of RS Pup, Westerlund derived a distance of 1 740 pc and 
an age of 4 million years for this stellar aggregate. Though this 
distance implies a spatial coincidence of RS Pup with Pup III, 
the age of the Cepheid clearly differs from that of the OB asso- 
ciation. 

The age of a Cepheid can be reliably estimated from the 
period-age relationship, whose existence is a consequence of the 
period-luminosity and mass-luminosity relationships, even if the 
actual crossing number of the Cepheid through the instability 
strip is uncertain. The up-to-date form of the period-age rela- 
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Table 5. Observationally determined distances of RS Pup. 
"IRSB" stands for "infrared surface brightness". 



Dist. [pc] Method Source 



1 780 ± 200 


Light echoes 


T T 1 i* 1 ^\T^ k 

Havlen l|1972a|l 


2040: 


Hipparcos parallax 


ESA (I1997II 


2iir^5 


IRSB 


Fouque et al. (|2003|l 


2 052^^^ 


IRSB 


Storm et al. (i2004| 


2 004 + 59 


IRSB (least squares) 


Barnes et al. (j2005l( 


1984+ 122 


IRSB (Bayesian) 


Barnes et al. (|2005|( 


CT 1 +266 


Hipparcos re-reduction 


van Leeuwen ([2007ab 


694+381 

"=^^-181 


Hipparcos re-reduction 


van Leeuwen et al. (12007b 


1 830+^f 


IRSB 


Fouque et al. (|2007|( 



1 992 + 28 Light echoes Present work 



tionship is (Bono et al. l2005l l: 

logf = (-0.67 + 0.01) logP-H (8.31 +0.08) 

where t is the age in years. This formula, valid for the average 
metallicity of our Galaxy, gives an age of about 17 million years 
for RS Pup excluding the possibility of physical membership of 
the Cepheid in a 4 million year old stellar grouping. The con- 
tradiction between the ages of the Cepheid and the association 
has already been noted by Havlen (1972b), though the age of 
the Cepheid was much underestimated using the period-age re- 
lationship available to him in the early 1970s. 

Later studies even cast doubt on the existence of the 
Pup III association (Eggen |l977l l. In the meantime Herbst d 19751 1 
pointed out the existence of an R association (i.e. an association 
embedded in reflection nebulosity) around RS Pup. This Pup R3 
association was found to be at a distance of 1 600 + 200 pc from 
us, possibly indicating a physical relation between the alleged 
Pup III and Pup R3 associations. It is problematic, however, that 
this R association was defined on only 3 arbitrarily selected stars 
whose distance moduli happened to be in accordance with that 
of RS Pup obtained by Havlen ( 1972al. 

The exhaustive list of the OB associations also includes 
PupOB3 but a distance of 2490pc is assigned to it (Mel'nik & 
Efremov 1995) . The source of this distance has not been given. 
In a more recent study, Bhatt et al. (1998) analyzed the polar- 
ization measurements in the direction of PupOB3. Their con- 
clusion was that the association members can be distinguished 
from the field stars based on their polarization characteristics. 
Bhatt et al. (1998) did not attempt to derive a new distance for 
PupOB3, instead they used Westerlund's d 19631 ) value of about 
1700pc. 

These contradictory statements and especially the age crite- 
rion lead us to be cautious, therefore we refrain from using the 
cluster/association method for deriving the distance of RS Pup. 

6. Conclusion 

We measured a geometric distance of li = 1 992 + 28 pc to 
the long-period Galactic Cepheid RS Pup using the propaga- 
tion of light echoes in its circumstellar nebula. Considering the 
intensity-mean magnitudes = 8.446 and my - 7.009 from 
Berdnikov (private communication), an extinction of E(B -V) - 
0.457 + 0.009 (Fouque et al. 2007 ), and the reddening laws from 
Laney & Stobie (119931 ), the absolute magnitudes of RS Pup in 
B and V are Mb = -5.11 + 0.05 and My = -6.09 + 0.05. Our 



distance is in good agreement with several other estimates, in- 
cluding the value derived by Havlen (1972a) of 1 780 + 200 pc 
using the same method. This is, to date, the most accurate ge- 
ometric distance to a Cepheid. A similarly straightforward dis- 
tance determination using both geometric and physical methods 
is possible in the far-infrared spectral region. The infrared echo 
is stimulated by the variable heating of the dust grains caused 
by the periodic temperature variations during the pulsation of 
RS Pup. The advantage of the infrared band is that the effect of 
the interstellar extinction can be circumvented and the method is 
independent of the scattering function of the dust grains. Though 
this method has been elaborated by Mayes et al. (119851) , it has 
never been applied until now. 

The light echoes of RS Pup are arguably the most spectac- 
ular example of this fascinating phenomenon. In contrast with 
supernovae or stellar outbursts (see e.g. Tylenda 2004), they 
are permanently observable using even moderately large tele- 
scopes, thanks to the large angular extension of the nebula (x; 
2.5 arcmin). In our images, we can monitor at all times the prop- 
agation of more than seven echoes, as visible on the interpolated 
movie available from the A&A web site. This nebula undoubt- 
edly holds a wealth of information about the mass-loss history 
of this star. It will thus be instrumental to understand the evolu- 
tion of Cepheids. Less than five years away from the centenary 
of the discovery of the Period-Luminosity relation by Leavitt & 
Pickering (119121 ), RS Pup could well become a "Rosetta stone" 
for this important class of stars. 
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Fig. 6. PSF subtracted images of RS Pup, sequentially in phase from the upper left to the lower right, 
for aU images. 



;. The intensity scale is identical 
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Fig. 8. Photometry of the ten nebular features used for the determination of the distance to RS Pup, with the corresponding adjusted 
photometric curves. The coordinates of each pixel in the images are indicated as "jc - y" in the upper left corner (See Table [3] for 
further details). The horizontal axis is in phase units, and the vertical axis is in analog-digital units (ADUs). 
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Fig. 1. Snapshots of a spline interpolation of our observations of RS Pup. An MPEG movie showing the progression of the wave- 
fronts in the nebula is available through the A&A web site. A ghost image of part of the telescope pupil can be seen on the right of 
the images. 



